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Abstract. Reeb graphs provide a method for studying the shape of a manifold by en- 
coding the evolution and arrangement of level sets of a simple Morse function defined on 
the manifold. Since their introduction in computer graphics they have been gaining pop- 
ularity as an effective tool for shape analysis and matching. In this context one question 
deserving attention is whether Reeb graphs are robust against function perturbations. Fo- 
cusing on 1 -dimensional manifolds, we define an editing distance between Reeb graphs of 
curves, in terms of the cost necessary to transform one graph into another. Our main result 
is that changes in Morse functions induce smaller changes in the editing distance between 
Reeb graphs of curves, implying stability of Reeb graphs under function perturbations. 



Introduction 

The shape similarity problem has since long been studied by the computer vision com- 
munity for dealing with shape classification and retrieval tasks. It is now attracting more 
and more attention also in the computer graphics community where recent improvements 
in object acquisition and construction of digital models are leading to an increasing ac- 
cumulation of models in large databases of shapes. Comparison of 2D images is often 
dealt with considering just the silhouette or contour curve of the studied object, encoding 
shape properties, such as curvature, in compact representations of shapes, namely, shape 
descriptors, for the comparison. The same approach is more and more used also in com- 
puter graphics where there has been a gradual shift of research interests from methods of 
representing shapes toward methods of describing shapes of 3D models. 

Since |24), Reeb graphs have been gaining popularity as an effective tool for shape anal- 
ysis and description tasks as a consequence of their ability to extract high-level features 
from 3D models. Reeb graphs were originally defined by Georges Reeb in 1946 as topo- 
logical constructs |22|. Given a manifold ./# and a generic enough real-valued function 
/ defined on the simplicial complex defined by Reeb, conventionally called the Reeb 
graph of (./#, /), is the quotient space defined by the equivalence relation that identifies the 
points of ^ belonging to the same connected component of level sets of /. Reeb graphs 
effectively code shapes, both from a topological and a geometrical perspective. While the 
topology is described by the connectivity of the graph, the geometry can be coded in a va- 
riety of different ways, according to the type of applications the Reeb graph is devised for, 
simply by changing the function /. Different choices of the function yield insights into the 
manifold from different perspectives. The compactness of the one-dimensional structure, 
the natural link between the function and the shape, and the possibility of adopting differ- 
ent functions for describing different aspects of shapes and imposing the desired invariance 
properties, have led to a great interest in the use of Reeb graphs for similarity evaluation. In 
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ifTSl . Hilaga et al. use Multiresolution Reeb Graphs based on the distribution of geodesic 
distance between two points as a search key for 3D objects, and the similarity measure 
constructed in this setting is found to be resistant to noise. In this approach resistance to 
changes caused by noise essentially relies on the choice of the geodesic distance to build 
the Reeb graph. In H, Biasotti et al. base the comparison of Extended Reeb Graphs on 
a relaxed version of the notion of best common subgraph. This approach gives a method 
for partial shape-matching able to recognize sub-parts of objects, and can be adapted to the 
context of applications since there is no requirement on the choice of the function /. Both 
ifTSll and PI present algorithms for similarity evaluation. 

To the best of our knowledge, mathematical assessment of stability against function 
perturbations is still an open issue as far as Reeb graphs are concerned. This question 
deserves attention since it is clear that any data acquisition is subject to perturbations, noise 
and approximation errors and, if Reeb graphs were not stable, then distinct computational 
investigations of the same object could produce completely different results. This paper 
aims to be possibly the first positive answer to this question. 

We confine ourselves to consider Reeb graphs of curves. In this setting Reeb graphs 
are simply cycle graphs with an even number of vertices corresponding alternatively to the 
maxima and minima of the function. We also equip vertices of Reeb graphs with the value 
taken by the function at the corresponding critical points. 

Our main contribution is the construction of a distance between Reeb graphs of curves 
such that changes in functions imply smaller changes in the distance. Our distance is based 
on an adaptation of the well-known notion of editing distance between graphs [25 1. We 
introduce three basic types of editing operations, represented in Table [H corresponding to 
the insertion (birth) of a new pair of adjacent points of maximum and minimum, the dele- 
tion (death) of such a pair, and the relabelling of the vertices. A cost is associated with 
each of these operations and our distance is given by the infimum of the costs necessary 
to transform a graph into another by using these editing operations. Our main result is the 
global stability of labelled Reeb graphs under function perturbations (Theorem l6.3b : 

Main Result. Let f,g:S^ — > M fee two simple Morse functions. Then the editing distance 
between the labelled Reeb graph of{S^, f) and that of{S^, g) is always smaller or equal to 
the C^-norm off — g. 

The main idea of the proof is to read editing operations in terms of degenerate strata 
crossings of the space of smooth functions stratified as in |6|. We also obtain a lower 
bound for our editing distance. Indeed, we find that it can be estimated from below by the 
natural pseudo-distance between closed curves studied in 1 13 1. 

The paper is organized as follows. In Section [T] we review some of the standard facts 
about Morse functions, the C topology, the theory of stratification of smooth real valued 
functions, and Reeb graphs. Section|2]deals with basic properties of labelled Reeb graphs 
of closed curves. Section [3] is devoted to the definition of the admissible deformations 
transforming a Reeb graph into another, the cost associated with each kind of deformation, 
and the definition of an editing distance in terms of this cost. Section |4] is intended to 
provide a suitable lower bound for our distance, the natural pseudo-distance; this represents 
a useful tool both to show the well-definiteness of our distance and to compute it in some 
simple cases. In Sections |5]and|6]it is shown that our distance is both locally and globally 
upper bounded by the difference, measured in the C^-norm, between the functions defined 
on 5^ Eventually, a brief discussion on the results obtained concludes the paper. 
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1. Preliminary notions 

In this section we recall some basic definitions and results about Morse functions and 
Reeb graphs. Moreover, with the aim of proving stability of Reeb graphs under function 
perturbations in mind, we recall some concepts concerning the space of smooth real valued 
functions on a smooth manifold: the topology and the theory of the natural stratification. 

Throughout the paper, ^ denotes a smooth (i.e. differentiable of class C°°) compact 
n-manifold without boundary, and ^(^,R) the set of smooth real functions on ^ . 



1.1. Simple Morse functions. Let us recall the following concepts from fT9l. 

Let / e ^(^,M). A point p G ^ is called a critical point of / if, choosing a local 
coordinate system (xi , . . . in a neighborhood U of p, it holds that 

and it is called a regular point, otherwise. Throughout the paper, we set — {p <^ ^ : 
p is a critical point of /}. 

If /? e K{f), then the real number f{p) is called a critical value of /, and the set {q e 
^ : q G / is called a critical level of /. Otherwise, if p ^ K{f), then f{p) is 

called a regular value. Moreover, a critical point p is called non-degenerate if and only if 
the second derivative matrix 

dxidxj 

is non-singular, i.e. its determinant is not zero. 

By the well-known Morse Lemma, in a neighborhood of a non-degenerate critical point 
p, it is possible to choose a local coordinate system (xi , . . . ,x„) such that 

.2 I „2 , I „2 



/-/(z') -4- •■■-4+4+1 + 



The number k is uniquely defined for each critical point p and is called the index of p. Such 
an index completely describes the behavior of / at p. For example, k means that the 
corresponding p is a minimum for f; k ~n means that /? is a maximum; < A: < n means 
that p is a saddle point for /. 

Definition 1.1. A function / G J^(^, K) is called a Morse function if all its critical points 
are non-degenerate. Moreover, a Morse function is said to be simple if each critical level 
contains exactly one critical point. 

It is well-known that every Morse function has only finitely many critical points (which 
are therefore certainly isolated points). The importance of non-degeneracy is that it is 
the common situation; indeed, in a sense that will be explained in Subsection 11.31 the 
occurrence of degenerate critical points is really quite rare. 

1.2. The C topology on the space of real valued functions. To topologize 
let us recall the definition of C'-norm, with < r < °o (see, e. g., ll20l I2TI '). Let {Ua} 
be a finite coordinate covering of with coordinate maps ha : Ua ^ , and consider 
a compact refinement {Ca} of {Ua} (i.e. Ca C Ua for each a, and [JCa = For 
/G ^(^,R), let us set/a = /o/i„i : ha{Ca)^K.. Then the C -norm of f is defined as 



ll/llc- =max<( max \fa{u)\, max 



dUj 



max 

u£ha(Ca 



je{\,...,n] Ji,...jre{l,. 



duj^ ■ ■ ■ duj^ 
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The above norm defines a topology on ^(^,M), known as the C topology (or weak 
topology), with < r < 00 (cf. 1 16] chap. 2]). In the following, we will denote by B,-{f,5), 
< r < 00, the open ball with center / and radius 5 in the C topology, i.e., g G Br{f, 5) 
if and only if ||/ — §11^/ < d. The C°° topology is simply the union of the topologies on 
.^{^, R) for every < r < 00. 

1.3. Natural stratification of tlie space of real valued functions. Let us endow ^(^, R) 
with the C°° topology, and consider the natural stratification of such a space, as exposed 
by Cerf in [Q (see also 12311 ). The natural stratification is defined as a sequence of sub- 
manifolds of ^ , R) , , ^ of co-dimension , 1 , . . . , y , . • • ^ respectively, 
that constitute a partition of , M), and such that the disjoint union U U . . . U 

is open for every j. 

Before providing a brief description of the strata, let us recall the following equivalence 
relation that can be defined on ^(^,R). 

Definition 1.2. Two functions /, g G R) are called topologically equivalent if there 

exists a diffeomorphism ^ : — > ^ and an orientation preserving diffeomorphism rj : 

M^Rsuchthat^((^(/7)) = r\{f{p)) for every p G 

The above relation is also known as isotopy in 16], and left-right equivalence in O. 

Let us describe and pointing out their main properties that allow us to leave 
aside the remaining strata. 

• The stratum is the set of simple Morse functions. 

• The stratum is the disjoint union of two sets and open in where 

- is the set of functions whose critical levels contain exactly one critical 
point, and the critical points are all non-degenerate, except exactly one. In 
a neighborhood of such a point, say p, a local coordinate system (xi, . . . ,x„) 
can be chosen such that 

/ = f{p) -x\-...-xl + . . . +xl. 

- is the set of Morse functions whose critical levels contain at most one 
critical point, except for one level containing exactly two critical points. 

is dense in the space ^(^,R) endowed with the C topology, 2 < r <°° (cf. 
lfT6l chap. 6, Thm. 1.2]). Therefore, any smooth function can be turned into a simple 
Morse function by arbitrarily small perturbations. Degenerate critical points can be split 
into several non-degenerate singularities, with all different critical values (Figure [T| (a)). 
Moreover, when more than one critical points occur at the same level, they can be moved 
to close but different levels (Figure[T] {b)). 

It is well-known that two simple Morse functions are topologically equivalent if and 
only if they belong to the same arcwise connected component (or co-cellule) of |6 , p. 
25]. 

is a sub-manifold of co-dimension 1 of J^*^ U^', and the complement of U^' 
in is of co-dimension greater than 1. Consequently, given two functions G we 
can always find f,g€ arbitrarily near to f,g, respectively, for which the path 

h{X) = {1-X)f + Xg, with A G [0, 1], is such that 

(1) and /, ^ are topologically equivalent to /, respectively; 

(2) /z(A) belongs to U for every A G [0, 1]; 

(3) /j(A) is transversal to J^'. 
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Figure 1 . («) a function / e admitting a degenerate critical point p (center) 
can be perturbed into a simple Morse function fi with two non-degenerate critical points 
p' ,p" (left), or into a simple Morse function /t without critical points around p (right); (b) 
a function / e (center) can be turned into two simple Morse functions /i ,/2, that are 
not topologically equivalent (left-right). 

As a consequence, h{X) belongs to for at most a finite collection of values A, and does 
not traverse strata of co-dimension greater than 1 (see, e.g., |[T4l ). 

1.4. Reeb graph of a manifold. In this subsection we restate the main results concerning 
Reeb graphs, starting from the following one shown by Reeb in ll22l . Here we consider 
pairs { J^ J), with connected and / G C ^(.-#,E). 

Theorem 1.3. The quotient space of ^ under the equivalence relation "p and q belong 
to the same connected component of the same level set of f" is a finite and connected 
simplicial complex of dimension 1. 

This simplicial complex, denoted by Yf, is called the Reeb graph associated with the 
pair (^,/). Its vertex set will be denoted by V(r/), and its edge set by E{r f). Moreover, 
if vi , V2 G V(r^) are adjacent vertices, i.e., connected by an edge, we will write e(vi , V2) G 
^(r^). Since the vertices of a Reeb graph correspond in a one to one manner to critical 
points of / on the manifold (see, e.g., 15, Lemma 2.1]), we will often identify each 
V G V{Tf) with the corresponding p G K{f). 

Given two topologically equivalent functions f,g€ it is well-known that the asso- 
ciated Reeb graphs, Ff and Fg, are isomorphic graphs, i.e., there exists an edge-preserving 
bijection <t> ; V(r/) ^ V{Fg). Beyond that, an even stronger result holds. Two functions 
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f,g G ^ are topologically equivalent if and only if such a bijection <I> also preserves the 
vertices order, i.e., for every v,w G V{rf), /(v) < f{w) if and only if g{<i>{v)) < g{<P{w)). 

The preceding result has been used by Arnold in |2 1 to classify simple Morse functions 
up to the topological equivalence relation. 

2. Labelled Reeb graphs of closed curves 

This paper focuses on Reeb graphs of closed curves. Hence, the manifold ^ that will 
be considered from now on is 5', and the function / will be taken in C ^(5\R). 
The Reeb graph Ff associated with is a cycle graph on an even number of vertices, 

corresponding, alternatively, to the minima and maxima of / on \2l \ (see, for example, 
Figure|2](fl) — (b)). Furthermore, we label the vertices ofFf, by equipping each of them 
with the value of / at the corresponding critical point. We denote such a labelled graph by 
(r/,/j ), where /j : V (F/) M is the restriction of / : 5' -> R to K{f). A simple example 
is displayed in Figure |2] (a) — (c). To facilitate the reader, in all figures of this paper we 
shall adopt the convention of representing / as the height function, so that (v^) < (v^) 
if and only if v„ is lower than v/, in the picture. 



/ 




(a) ib) (c) 

Figure 2. (a) a pair (S',/), with /the height function; (&) the Reeb graph associ- 
ated with (S',/); (c) the labelled Reeb graph (Fy,/j) associated with (5',/). Here labels 
are represented by the heights of the vertices. 

The natural definition of isomorphism between labelled Reeb graphs is the following 
one. 

Definition 2.1. We shall say that two labelled Reeb graphs (r^,/j ), {Fg,g^ ) are isomorphic 
if there exists an edge-preserving bijection <I> : V(r^) — V{rg) such that (v) = g^ (^>(v)) 
for every v G V(r/). 

The following Proposition 12.41 provides a necessary and sufficient condition in order 
that two labelled Reeb graphs are isomorphic. It is based on the next definition of re- 
parameterization equivalent functions. 

Definition 2.2. Let Jif{S^) be the set of homeomorphisms on SK We shall say that two 
functions f,g G C J^(5'',R) are re-parameterization equivalent if there exists T G 
J^{S^) such that /(/?) = g{T:{p)) for every p(GS\ 

Lemma 2.3. Let (Tf,/^) and (Tgjg^ ) be labelled Reeb graphs associated with {S^ ,f) and 
{S^ ,g), respectively. If an edge-preserving bijection <P : V(r/) — > V{rg) exists, then there 
also exists a piecewise linear T G J^{S^) such that T|^j^ ^ = «t>. If moreover — g^o^, 
then f = go T. 
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Proof. The proof of the first statement is inspired by 1.13. Lemma 4.2]. Let us construct 
T by extending <I> to 5' as follows. Let us recall that V{Tf) = K{f) and ViTg) =^ K{g), 
and, by abuse of notation, for every pair of adjacent vertices p' ,p" £ V{rf), let us iden- 
tify the edge e{p' ,p") G E{Tf) with the arc of 5^ having endpoints p' and p" , and not 
containing any other critical point of /. For every p G K{f ), let T{p) = ^(p). Now, 
let us define T{p) for every p G 5' \K{ f ). Given p G \ K{ f), we observe that there 
always exist p' ,p" G V{rf) such that p G e{p',p"). Since <I> is edge-preserving, there ex- 
ists e(<I> (/?'), <!>(/?")) — e{z{p'), i{p")) G E{Tg). Hence, we can define 'i{p) as the unique 
point of e(T(p'),T(/?")) such that, if f{p) = {\ ~ Xp)f{p')+Xpf{p"), with A,, G [0,1], then 
g{l[p)) = (1 — Xp)g{x{p')) + Xpg{T{p")). Clearly, T belongs to ^(5') and is piecewise 
linear. 

As for the second statement, it is sufficient to observe that, if /j = ^| o<I>, since T{p) = 
^{p) for every p G K[f), then clearly (p) — g^ {t{p)) for every p G K{f). Moreover, for 
every p G 5' \K{f), by the construction of T, it holds that gixip)) = (1 - Xp)g{(i>{p')) + 
^pgW)) = (1 - + ^pfip") = fip)- In conclusion, f{p) = g{x{p)) for every 

p G 5', and, hence, f,g are re-parameterization equivalent. 

□ 

Proposition 2.4 (Uniqueness theorem). Let (Ff,/^), {Fgjg^) be labelled Reeb graphs as- 
sociated with {S^,f) and {S^,g), respectively. Then (Tf,/^) is isomorphic to {Tgjg^) if and 
only iff and g are re-parameterization equivalent. 

Proof. The direct statement is a trivial consequence of Lemma 123] 

As for the converse statement, it is sufficient to observe that any T G ^(5') such that 
f — gox, as well as its inverse T~\ takes the minima of / to the minima of g and the 
maxima of / to the maxima of g. Hence, <I> : V(r/) V{Tg), with <I> — t^^^^ j, is an edge 

preserving bijection such that /j = g| o <I>. □ 

As a consequence of Proposition l2.4l two labelled Reeb graphs isomorphic in the sense 
of Definition l2.1l will always be identified, and in such case we will simply write (r/,/j) = 

The following Proposition 12.51 ensures that, for every cycle graph with an appropriate 
vertices labelling, there exists a unique (up to re-parameterization) pair (5^/), with / G 
having such a graph as the associated labelled Reeb graph. 

Proposition 2.5 (ReaUzation theorem). Let {G,£) be a labelled graph, where G is a cycle 
graph on an even number of vertices, and i : V{G) — >■ R is an injective function such that, 
far any vertex v^ adjacent {that is connected by an edge) to the vertices v\ and V3, either 
both £{vi) and £(v^) are smaller than i{v2), or both £{v\) and £{vj,) are greater than £(v2). 
Then there exists a simple Morse function / : 5' — > M such that i^fif^^yi ) ~ {Cr,£). 

Proof. It is evident. □ 

By virtue of the above Uniqueness and Realization theorems (Propositions |2.4| and |2.5l ), 
for conciseness, when a labelled Reeb graph will be introduced in the sequel, the associated 
pair will be often omitted. 

3. Editing distance between labelled Reeb graphs 

We now define the editing deformations admissible to transform a labelled Reeb graph 
of a closed curve into another We introduce at first elementary deformations and then the 
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deformations obtained by their composition. Next, we associate a cost with each type of 
deformation, and define a distance between labelled Reeb graphs in terms of such a cost. 

Definition 3.1. Let (r/,/j) be a labelled Reeb graph with 2n vertices, n > 1. We call an 
elementary deformation of {Tf^f^) any of the following transformations: 

(B) (Birth): Assume e(vi,V2) G E{Yf) with /j(vi) < /j(v2)- Then {TfJ^) is trans- 
formed into a labelled graph (G, (.) according to the following rule: G is the new 
graph on 2n + 2 vertices, obtained deleting the edge e(vi , V2) ™d inserting two new 
vertices mi, M2 and the edges e{vi,ui),e{ui,U2),e{u2,V2)', moreover, £ : V{G) — > K 
is defined by extending /j from V(r^) to V{G) — V(r^) U {ui,U2} in such a way 
that^H,(r^.) =/j,and/j(vi) <^(m2) <^(mi) </j(v2). 

(D) (Death): Assume « > 2, and e{vi,ui),e{ui,U2),e{u2,V2) G E{rf), with /j(vi) < 
/j ("2) < /j (mi) < /j (^2)- Then (r/,/j ) is transformed into a labelled graph {G,£) 
according to the following rule: G is the new graph on 2n — 2 vertices, obtained 
deleting ui, U2 and the edges e{vi,ui), e{ui,U2), e{u2,V2), and inserting an edge 
e(vi,V2); moreover, £ : V{G) R is defined as the restriction of to V(r^) \ 

{u\,U2}. 

(R) (Relabelling): (r/,/j ) is transformed into a labelled graph {G,£) according to the 
following rule: G = Tf, and for any vertex V2 adjacent to the vertices vi and V3 
(possibly vi = V3 for n — I), if both /j (vi) and (V3) are smaller (greater, respec- 
tively) than f^{v2), then both £{vi) and £{vj,) are smaller (greater, respectively) 
than £{v2)', moreover, for every v ^ w, £{v) ^ £(yi). 
We shall denote by T{rf,f^) the result of the elementary deformation T applied to (r^,/j ). 

Table [T] schematically illustrates the elementary deformations described in Definition 

mi 

Proposition 3.2. Let T be an elementary deformation of {Tf,f^), andlet (G,£) — T{Yf,f^). 
Then {G,£) is a Reeb graph (Fg^g^) associated with a pair {S^ ,g), and g £ is unique 
up to re-parameterization equivalence. 

Proof. The claim follows from Propositions 12.51 and 12 .41 □ 
As a consequence of the above result, from now on, we will directly write T{Tf,f^^ ~ 

Moreover, since the previous Proposition 13.21 shows that an elementary deformation 
of a labelled Reeb graph is still a labelled Reeb graph, we can also apply elementary 
deformations iteratively. This fact is used in the next Definition 13. 3 1 

Given an elementary deformation T of (r/^, /j) and an elementary deformation S of 
T{rf, f^ ), the juxtaposition ST means applying first T and then S. 

Definition 3.3. We shall call deformation of {rf,f^) any finite ordered sequence T = 
{Ti,T2, . . . ,Tr) of elementary deformations such that Ti is an elementary deformation of 
(Tfjf^ ), T2 is an elementary deformation of Ti (r/,/j ), Tr is an elementary deformation 
of Tr-iTr-2- ■ ■Ti{rf,f^). We shall denote by r(r/,/j) the result of the deformation T 
applied to (r^,/j). 

Let us define the cost of a deformation. 

Definition 3.4. Let T be an elementary deformation transforming {rf,f) into {rg,g.). 
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Table 1 . The upper two figures schematically show the elementary deformations of 
type (B) and (D), respectively; the third figure shows an example of elementary deforma- 
tion of type (R). 



• If r is of type (B) inserting the vertices Mi , M2 G V{rg), then we define the associ- 
ated cost as 



C{T): 



I^|(«i)-^|(m2)|_ 



2 

• If r is of type (D) deleting the vertices mi,M2 S V{rf), then we define the associ- 
ated cost as 

|/(«l)-/(«2)| 



c{T)-' ' 



2 

• If r is of type (R) relabelUng the vertices v G V{rf) = V'(rg), then we define the 
associated cost as 



c{T)= max |/j(v)-g|(v)|. 



10 



B. DI FABIO AND C. LANDI 



Moreover, if T = (Ti , . . . , Tr) is a deformation such that Tr - ■ ■Ti{rf,f^) = (Fgjg^), we 

r 

define the associated cost as c{T) — Y, c(7/)- 

(=1 

We now introduce the concept of inverse deformation. 

Definition 3.5. Let T be a deformation such that r(r/,/j ) = (Fgjg^). Then we denote by 
r^', and call it the inverse of T, the deformation such that T^^ (Tgjg^) — (Tf,/^ ) defined 
as follows: 

• If r is elementary of type (B) inserting two vertices, then T^^ is of type (D) 
deleting the same vertices; 

• If r is elementary of type (D) deleting two vertices, then T^' is of type (B) insert- 
ing the same vertices, with the same labels; 

• If r is elementary of type (R) relabelling vertices of y(r^), then T^^ is again of 
type (R) relabelling these vertices in the inverse way; 

. If r = (ri , . . . , ?;), then r-' = (r- rf' ). 

Proposition 3.6. For every deformation T such that T(rf,f^) — (Tgjg^), c(r^') = c[T). 

Proof. Trivial. □ 

We prove that, for every two labelled Reeb graphs, a finite number of elementary de- 
formations always allows us to transform any of them into the other one. We recall that 
we identify labelled Reeb graphs that are isomorphic according to Definition l2.1l We first 
need a lemma, stating that in any labelled Reeb graph with at least four vertices we can 
find two adjacent vertices that can be deleted. 

Lemma 3.7. Let {Tf,f^) be a labelled Reeb graph with at least four vertices. Then there 
exist e{v\,u\),e{u\,U2),e{u2,V2) eE{Tf), with f^{vi) </j(M2) < f^{ui) </j(i'2)- 

Proof. Let V(r^) — {ao,bo,ai,bi, . . .a,„_i,b,„_i}, m > 2. In the following, we convene 
that, for k ^Z, at: and b^ are equal to a^i^ modm) and ^(A- mod;?i)^ respectively. We assume 
that £'(r/') = {e{ai,bi) : ; > 0} U {e(^;,a,+i) : / > 0}, and /j (a,) < /j (fe,) for every From 
the definition of labelled Reeb graph associated with a pair (5' , /), it follows that /j (Z?,) > 

fl (ai+i), /] (a;) ^ /j (flr+i), /] (hi) ^ /] (bi+i), for every /. 

The claim can be restated saying that there is at least one index / such that either (/) 
</j(a,+i) and/j(/7,) </j(/7,+ i) or (//) (fl,+i) < /j (a,) and (^7,) < /j (Vi) hold. 
We prove this statement by contradiction, assuming that for every / > neither (/) nor 
(//) hold. Since (/) does not hold, either /j (oo) > /j (fli) or (bo) > f^{b\) or both. Let 
us consider the case when /j [bo) > (bi). Since (//) does not hold either, it follows that 
/j(fl2) >./j(«i)- Recalling that (/) does not hold, we obtain /j (fe i ) > /j (^2 ) ■ Iterating the 
same argument, we deduce that /j (/?,) > f^{bi + I) for every i > 0, contradicting the fact 
that bi„ — bo. An analogous proof works when we consider the case /j (aq) > (oi ). □ 

Proposition 3.8. Let (ry,/j ) and {Tg,g^ ) be two labelled Reeb graphs. Then the set of all 
the deformations T such that TiYf^f^) — {Yg^g^) is non-empty. This set of deformations 
will be denoted by ^{{Tf,f^),{Tg,g^)). 

Proof. If (F/j /j) — (Fg,g|), then it is sufficient to take the elementary deformation T of 
type (R) transforming (F ^,/j) into itself. Otherwise, if (Fj,/j ) ^ (ri,,g| ) and F^ has at least 
four vertices, by Lemma [STtI we can apply a finite sequence of elementary deformations 
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of type (D) to (Ff,/^), so that in the resuhing labelled Reeb graph (r/,,/!|), F/, has only two 
vertices, say m, v, with (m) < (y). If (rg,^|) has also at least four vertices, by Lemma 
13.71 there exists a finite sequence of elementary deformations of type (D) to (rs,,^|), say 
S — {Si,. . . ,Sp), so that in the resulting labelled Reeb graph (r/,/,/!|), F/,/ has only two 
vertices, say m',v', with /i'(m') < h'{v'). So, we can apply to (r/,,/Z|) an elementary defor- 
mation of type (R) so to obtain (F/,/ ,h'^. Finally, by Definition l3.5l we can apply to (F/,/ , h'^) 

the finite sequence of elementary inverse deformations of type (B), S^^ = {Sp\. ..,S^^), 
in order to obtain (Fg,g|). For (F/^,/j) or (Fg,§|) with only two vertices, the same proof 
applies without need of deformations of type (D) or (B), respectively. □ 

A simple example explaining the above proof is given in Figure[3] 




M4 


M4 




M6 


(B^\ 












Ml 


Ml 



Ml 



Figure 3. The leftmost labelled Reeb graph is transformed into the rightmost one 
applying first three elementary deformations of type (D), then one elementary deformation 
of type (R), and finally two elementary deformations of type (B). 



We point out that the deformation constructed in the proof of Proposition 13.81 is not 
necessarily the cheapest one, as can be seen in Example |2] 

We now introduce an editing distance between labelled Reeb graphs, in terms of the 
cost necessary to transform one graph into another. 

Theorem 3.9. Forevery two labelled Reeb graphs (Fy,/j) and {Tg,g^), we set 



.„,5i))= inf c{T). 



Then d is a distance. 



The proof of the above theorem will be postponed to the end of the following section. 
Indeed, even if the properties of symmetry and triangular inequality can be easily verified, 
the property of the positive definiteness of d is not straightforward because the set of all 
possible deformations transforming {Tf,f^) to (Fi,,g|) is not finite. In order to prove the 
positive definiteness of d, we will need a further result concerning the connection between 
the editing distance between two labelled Reeb graphs, (F/,/j), (Fg,^|), and the natural 
pseudo-distance between the associated pairs {S^,f), {S\g). 

4. A LOWER BOUND FOR THE EDITING DISTANCE 

Now we provide a suitable lower bound for our editing distance by means of the natural 
pseudo-distance. 

The natural pseudo-distance is a measure of the dissimilarity between two pairs {X, (p), 
{Y, y/), with X and Y compact, homeomorphic topological spaces and (p :X —^R, y/ .Y 
R continuous functions. Roughly speaking, it is defined as the infimum of the variation 
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of the values of (p and when we move from X to Y through homeomorphisms (see 
linilllKn] for more details). 

Such a lower bound is useful for achieving two different results. The first result, as 
mentioned in the preceding section, concerns the proof of Theorem 13.91 i.e., that li is a 
distance (see Corollary I4.2l i. The second one is related to an immediate question that 
can arise looking at the definition of d: Is it always possible to effectively compute the 
cheapest deformation transforming a labelled Reeb graph into another, since the number 
of such deformations is not finite? By using the natural pseudo-distance, we can estimate 
from below the value of d, and, in certain simple cases, knowing the value of the natural 
pseudo-distance allows us to determine the value of d (see, e.g., ExamplesfT}©. 

The following Theorem 14.11 states that the natural pseudo-distance computed between 
the pairs {S^ ,f) and (5' , is a lower bound for the editing distance between the associated 
labelled Reeb graphs. 

Theorem 4.1. Let (Ff,/^), (r^j^i) be labelled Reeb graphs associated with (5',/) and 

{S\g), respectively. Thend{{TfJ),{Tg,g )) > inf \\f-gox\\(^. 

' ' Te,jsr(si) 

Proof. Letusprovethat,foreveryre3^((rf,/;),(r„,gi)),c(r)> inf ||/-goT|U. 

First of all, assume that T is an elementary deformation transforming {Tf,f^^ into 
(Egj^l). For conciseness, slightly abusing notations, we will identify arcs of 5^ having 
as endpoints two critical points p',/?" e V (F/), and not containing other critical points of 
/, with the edges e{p',p") G E{Yf). 

(1) Let T be of type (R) relabelling vertices of y(F/). Since, by Definition 13. II (R), 

Tf~ Tg, we can always apply Lemma 123] considering <I> as the identity map, to 

obtain a piecewise linear T G J>^{S^) such that x{p) — p for every p e K{f). As 

far as non-critical points are concerned, following the proof of Lemma 12.31 for 

every p \K{f), x{p) is defined as that point on 5' such that, if /5 G e{p\p") G 

E{Tf), with fip) - (1 - Xp).np')+Xpf{p"), Xp G [0, 1], then T{p) G e{p\p") 

with g{T{p)) = (1 — Xp)g{p') + Xpg{p"). Therefore, by substituting to f{p) and 

g(T(/7)) the above expressions, we see that max|/(/7) —g(T(/7)) I = max \f,(p) — 

pes' peviFf) 

g,{p)\^c{T). 

(2) Let r be of type (D) deleting ^ i , ^2 e ^ (F/ ), the edges e(p\,qi),e{q\,q2),e{q2,P2), 
and inserting the edge e{p\,p2). Thus, for every p G A'(/)\{^i,^2}, f{p) — 
g{p). It is not restrictive to assume that f{p\) < fiqi) < f{qi) < fipi)- Then 
we can define a sequence (t„) of piecewise linear homeomorphisms on 5' ap- 
proximating this elementary deformation. Let Z„{p) = p for every every p G 
V{rf)\{qi,q2} — V(Fg) and n G N. Moreover, let ^be the point of e(/? 1,752) G 
E{Tg) such that g(^) = (such a point ^ exists because ^(/^i) — f{pi) < 
fill) < f{q\) < fipi) = gipi) and it is unique because we are assuming that 
no critical points of g occur in the considered arc). Let us fix a positive real 
number c < mm{g{p2) —g{q),giq) — gipi)}- For every n G N, let us define 
T^nili) (resp. Tniqi)) the only point on belonging to the arc with endpoints 
pi,q (resp. q,p2) contained in e{pi,p2), such that g{Tn{q\)) = g{q) - (resp. 
gi'^niqi)) = §® + f ) as shown in Figure H) Now, let us linearly extend T„ to 
all 5' in the following way. For every p € K{f), if p belongs to the arc 
with endpoints p' ,p" G K{f) not containing any other critical point, and is such 
that /(/?) = (1 ~Xp)f[p')+Xpf{p"), Xp G [0, 1], then x„{p) belongs to the arc 
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with endpoints x„[p')^ "^nip") not containing any other critical point, and is such 

that g{l„{p)) = (1 - Xp)g{t„{p')) +Xpg{T„{p")). Hence, T„ is piecewise linear 

for every n e N, and limmax|/(/9) — g(T„(;?))| — lim max \f{p) ~ g{'^n{p))\ ~ 
"^'^pes^ n^°°pev{rf) 

lim max{/(^i) -g(T„(^i)),/(^2) -g(T„fe))} = \f{qi)-g{q)\ = 5^ = 

c{T). 




V 



q2 

J J 

P\ '^n{P\)=P\ 

Figure 4. The construction of the homomorphism T„ as described in step (2) of the 
proof of Theorem l4.1l The arc e{p\,qi ) (e(q\,q2), and e{q2,P2), respectively) is piecewise 
linearly taken to the ai'c having T„(pi),T„((/i) (T„(gi), T„(g2) and T„(^2),T„(p2). respec- 
tively) as endpoints. 



(3) Let T be of type (B) deleting e{pi,p2) G E{rf), and inserting two vertices qi,q2 
and the edges e{pi,qi), e{q\,q2), e{q2,P2)- Then we can apply the same proof as 
(2), by considering the inverse deformation T^' that, by Definition [33] is of type 
(D) and, by Proposition 13. 61 has the same cost of T. 



Therefore, observing that in (1), the piecewise linear T can be clearly replaced by a se- 
quence (t„), with T„ = T for every n G N, we can assert that, for every elementary defor- 
mation T, there exists a sequence of piecewise linear homeomorphisms on S^, (t„), such 
thatc(r) = lim||/-goT„||co > inf ||/-^ot||co. 

Now, let T = {Ti,...,Tr) G £^{{rf,f^),{rg,g^)) and prove that, also in this case, c{T) > 
inf \\f-goT\\^. Let us set Tr--Ti{rf,f ) = (r^(,),/»), / = /W, g = For 

i— 1 , . . . , r, let (ri'^ be a sequence of piecewise linear homeomorphisms on for which 
it holds that c{Ti) = lim -/<') o II (--0, and let (Ih"'')^ be the constant sequence 
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such that ri"^ = Id for every n G N. Then 



c(r) = i:c(7^0 = lim||/«)-/''oTi')||co + EHm||/»-/('+%Tr 



X lim 11/^'' o ri'' o . . . o - /<'+') o 7,1'+') o ^ o-.-o | 



> lim||/W-/WoTi'''oTi'-''o...oTi'')||co> inf ||/-goT|lco, 
where the third equality is obtained by observing that 

/(') o T« o . . . o rf -/('■+') o Ti'+" o ^ o . . . o = (/(') o Ti'+^') o ri'^ o . . . o rf 

for every / G { 1 , . . . , r — 1 }, and that || • 11^0 is invariant under re-parameterization; the first 
inequahty is consequent to the triangular inequaUty. □ 

CoroUary4.2. Ifd{{TfJ^),{Tg,g^))^Qthen (F^,^) = (r,,^,). 

Proof. From Theoremgll J((r/,/),(r„,g|)) =0 implies that inf \\f ~ goiW^) =Q. 

In i?) it has been proved that when inf 1 1 f — ^ o t 1 1 ro = 0, with X, Y two closed curves 

te,^{x.Y) ^ 

of class at least C^, a homeomorphism T G J^{X , Y) exists such that f — goT. Therefore, 
the claim follows from Proposition 12. 41 □ 



Proof of Theorem \3.9\ The positive definiteness of d has been proved in Corollarv 14.21 
the symmetry is a consequence of Proposition [32J the triangular inequality can be easily 
verified in the standard way. □ 

Now we describe two simple examples showing how it is possible to compute the editing 
distance between two labelled Reeb graphs, (Fy^, /j ), (Fg,g| ), by exploiting the knowledge 

of the natural pseudo-distance value between the associated pairs (5^ /) , (5' , ^) . In partic- 
ular. Example [Uprovides a situation in which the infimum cost over all the deformations 
transforming (Fy^,/j) into {Tg,g^) is actually a minimum. In Example |2] this infimum is 
obtained by applying a passage to the limit. 

Example 1. Let us consider the two pairs (5' ,/), (5' depicted in Figure|5l with f,g€ 
We now show that ii((F/,/j), (F?,g|)) = ^if{qi) - fipi))- Indeed, in this case, the 
natural pseudo-distance between (5^,/) and {S^ ,g) is equal to j{f{qi) ^ f{pi)) (cf. |13]). 
Therefore, by TheoremgU it follows that t/((F/,/j ), (Fg,g|)) > \{f{qi)-f{pi))- On the 
other hand, the deformation T of type (D) that deletes the vertices pi,qi G V (F/), the edges 
e{p,qi),e{qi, pi),e{pi,q) and inserts the edge e{p,q) transforms (Fy, /j ) into (Fg,g| ) with 

cost c(r) = ^{f{qi)~f{pi)). Hence ^((F^,^), (Fg,^,)) = ^{f{qi)~f{pi)). 

Example 2. Let us consider now the two pairs (5' ,/), {S^ ,g) illustrated in Figure|6l Let 

fiqi) ^ f{P\) = /(^2) — fipi) = a. Then, clearly, inf ||/~go T||ro — %. Let us show 

that the editing distance between {Tf,f^) and (F|,,^|) is |, too. For every < £ < |, we 
can apply to (F/,/j) a deformation of type (R), that relabels pi,p2,qi,q2 in such a way 
that f{pi) is increased of | — e, and f{qi) is decreased of | — e for ; = 1,2, composed 
with two deformations of type (D) that delete pi with qi, i — 1,2. Thus, since the total 
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Figure 5. The pairs considered in Example [T] In this case rf({ry,/j ), (rg,g| )) = 
inf ||/-goT||co = |(/(9i)-/(/>i)) 

{S\f) q q' {S\g) 




Figure 6. The pairs considered in Example[2] Even in this case d[{T f ,f^),{Tg,g^)) = 
inf \\f-goz\\^ = {(f(q,)-f{pi)) 



cost is equal to | — e + 2e, by the arbitrariness of e, it holds that t/((r/,/j), (rg,g| )) < |. 
Applying Theorem |4.1[ we deduce that ii((r^,/j), (rg,^|)) = |. 

5. Local stability 

This section is intended to show that labelled Reeb graphs of closed curves are stable 
under small function perturbations with respect to our editing distance (see Theorem |5.5l ). 
The main tool we will use is provided by Theorem 15. 31 that ensures the stability of sim- 
ple Morse function critical values. This latter result can be deduced by the homological 
properties of the lower level sets of a simple Morse function / on a manifold and 
its validity does not depend on the dimension of Therefore, it will be given for any 
smooth compact manifold without boundary. 

For every / e ^(^,R), and for every a G M, let us denote by /" the lower level set 
(— oo, a] = {p E ^ : f{p) < a}. Let us recall the existing link between the topology of 
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a pair of lower level sets (/'',/"), with a,b & a < b, regular values of /, and the critical 
points of / lying between a and b. The following statements hold (cf. 1191 ): 

(St. 1) If the interval f^^{[a,b]) contains no critical points, then /" is a deformation re- 
tract of f\ so that the inclusion map /" ^ is a homotopy equivalence. 

(St. 2) If f^^{[a,b]) contains exactly one critical point of index k, then, denoting by G 
the homology coefficient group, it holds that 



G, ifk^k 
0, otherwise. 



In the remainder of this section we require / to be a simple Morse function. Accordingly, 
it makes sense to use the terminology critical value of index k to indicate a critical value 
that is the image of a critical point of index k. 

Lemma 5.1. Let f £ C ^(.^,M), and let a,b eR, a < b, be regular values off. If 
there exists A; G Z such that Hj{f^ ,f") ^ 0, then [a,b] contains at least one critical value 
of index k. 

Proof. From (St. 1), the absence of critical values in [a,b] implies that the homomorphism 
induced by inclusion i^. : Hk{f") Hkif'') is an isomorphism for each A; e Z. Conse- 
quently, by using the long exact sequence of the pair: 

— > H.if) ^ A H.ifj") -% H,_,{f") n ^ ... , 

it is easily seen that, for every k ^"L, the surjectivity of /<. and the injectivity of imply 
the triviality of Hk{f'\f"). This proves that if there exists A^ e Z such that Hj(^ f'',f") ^ 
0, then [a,b] contains at least one critical value of /. That the index of at least one of 
the critical values of / contained in [a,b] is exactly k is consequent to the sub-additivity 
property of the rank of the relative homology groups and to (St. 2). In fact, let ci , . . . , c„, be 
the critical values of / belonging to [a,b], and let so,.- - ,Sm be m + 1 regular values such 
thata = io < ci < ii < C2 < . . . <i,„-i <c,„ <5,„ ^b. Since it holds that rank//j(/'',/") < 

m 

^ rankf/^(/*' ), and by hypothesis rank//^(/ ,/") > 1, there exists at least one index 

(=1 

/ e { 1 , . . . , m} such that Hj{f''' , f'- ' ) 7^ 0. Now, applying (St. 2) with a replaced by 1 
and b replaced by s,-, we deduce that c, is a critical value of / of index k. □ 

The above statements (St. 1-2), Lemma |5T1 together with the following lemma, that is 
a reformulation of Lemma 4. 1 in ifTTl , provide the tools for proving the stability of critical 
values under small function perturbations (Theoi'em |5.3t . 

Lemma 5.2. Let Xi,X2,Xi, X[ , ^2 ' ^3 topological spaces such that Xi C X2 C C X[ C 
^2 ^ ^3. Hki^3,^i) = 0, Hf;{X^,Xi) = for every k ^1,. Then the homomorphism 
induced by inclusion Hi^{X[,Xi) — > Hk{X2-,X2) is injective for every A: G Z. 

Theorem 5.3 (Stability of critical values). Let f e C ^(.-#,R) and let c be a critical 
value of index k of f. Then there exists a real number 8{f,c) > such that each g G 
verifying \\f ~ g\\(0 < 5, < 5 < 5{f,c), admits at least one critical value of index k in 
[c~-5,c + 5]. 

Proof. Since / is Morse, we can choose a real number 5{f,c) > such that [c — 3 ■ 
5{f,c),c + 3 ■ 5(/,c)] does not contain any critical value of / besides c. Let < 5 < 
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5(/,c), and let g be a simple Morse function such that ||/ — ^||(-.o < 5. If 5 = 0, then the 
claim immediately follows. Let 5 > 0. Then, for every n e N, 

Since [c — 5 ■ , c — 5 /n] and [c + 5/n,c + 5 ■ ^^iti ] do not contain any critical value 

of / for every « e N, both Hk{f-^/" and //a:(/'+^'^ J'+^Z") are trivial for 

every A; G Z, and « G N. Consequently, from Lemma 15.21 the homomorphism induced 
by inclusion Hk{f+^/",f-^-^) Hk{g^+^-'^ ,g^-^-'^) is injective for each k e Z, 
and n G N. Moreover, since, for every n G N, [c — 5 • ^iiti + d/n] contains c, that is 

a critical value of index k of /, from (St. 2), it holds that Hj{f+^l'\f-^-'"^) ^ for 
every n G N. This fact, together with the injectivity of the above map, implies that also 
H-j^{g'^'^^'^ Tg'^^^'^) ^ for every n G N. So, by Lemma ISTl for every n G N, there 
exists at least one critical value c,', of index kof g with cj, G (c — 5 • ^^,c + 5 • ^^). By 
contradiction, let us suppose that [c — 5,c + 5] contains no critical values of index k of 
g. Then, since g is Morse, there would exist a sufficiently small real number e > such 
that (c — 5 — e,c + 5 + e) does not contain critical values of index kof g either, giving an 
absurd. □ 

We now prove the local stability of labelled Reeb graphs of closed curves. We need a 
lemma that holds for manifolds of arbitrary dimension. The global stability will be exposed 
in the next section. 

Lemma 5.4. Let f G C ^(^,M). Then there exists a positive real number 5{f) 
such that, for every 5, < 5 < 5{f), and for every g G r^", with \\f ~ g\\(j2 < 5, an edge 
and vertices order preserving bijection <I> : V{Tf) — > V(Tg) exists for which max \f (v) — 

g^mv))\<5. 

Proof. Let /:> i , . . . , p„ be the critical points of /, and ci,...,c„ the respective critical values, 
with c, < c,+ i for ; = l...,« — 1. Since is open in ^(^,M), endowed with the 

topology, there always exists a sufficiently small 5 (/) > 0, such that the closed ball 
with center / and radius 5(/), B2{f,5{f)), is contained in Moreover, 5{f) can be 
chosen so small that, for every ; = 1,...,« — 1, the intervals [c,- — 5(/),c,- + 5{f)] and 
[cj+i - 5(/),c,+i + 5{f)] ai-e disjoint. 

Fixed such a for every real number 5, with < 5 < 8{f), and for every g G 
such that II/ — g||c2 ^ 5, / and g belong to the same arcwise connected component of 
endowed with the C°° topology, and, therefore, are topologically equivalent functions. 
Consequently, there exists an edge and vertices order preserving bijection <I> : V{rf) 

y (r„) (see Subsection |L4t . Let us prove that <I> is such that max \f (v) — g. (<I>(v)) | < 5. 

veviTf) I 

Since / and g are topologically equivalent, it follows that g has exactly n critical points, 
p[,.. .,p'„. Let c'l = g{p[), . . . ,4 = g{p'„). We can assume < cj^j, for 1, ...,«- L 
The assumption ||/ — §||c2 < ^ implies that ||/ — g||co < 5. Therefore, by the previous 
Theorem 15.31 for every critical value c, of /, there exists at least one critical value of g 
of the same index of c, belonging to [c,- — 5,c,- + 5]. Moreover, since [c, — 5,c, + 5] fl 
[ci+i — 5,c,+ i + 5] — a for every i — 1, ...,« — 1, it follows that G [c,- — 5,c,- + 8] for 
every i — !,...,«. Hence, since <I> preserves the order of the vertices, necessarily 4>(/7,) ~ 

pj, yielding that max |/ (v) (<I>(v))| = max \f{pi) - g{<i>{pi))\ ^ ma\\ci ~ c-\ < 

veViTf) I I neKif) ' ' l<i<n 

5. a 
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Theorem 5.5 (Local stability). Let f G C ^(5 ,M). Then there exists a positive real 
number 5{f) such that, for every 5,Q<5 < 8{f), and for every g ^ .^^\ with ||/ — ^11^2 < 
5, it holds that d({T f,f^), {Tg,g^)) < 8. 

Proof. By Lemma |54l an edge and vertices order preserving bijection <I> : V{rf) — > V{rg) 

exists for which max \f{v) — g.{<P{v))\ < 5. Necessarily <I> takes minima into minima 

vev{rf) ' 

and maxima into maxima. Therefore, (Tfjg^ o't') = T{Tf,f^), with T an elementary de- 
formation of type (R), relabelling vertices of ViTA, having cost c{T) = max |/(v) — 

v&{Yf) 

gj(<I>(v))| < 5. Moreover, let us observe that ij^f.g^ o't') is isomorphic to (rg,g|) as la- 
belled Reeb graph (see Definition 12.1b . Thus, d{{rf,f^),{rg,g^)) ~ d{{r f ,f^)^{T f,g^ o 

4>)) = inf c(T) <5. □ 

re^{{r^,/|),(r„g|)) 

6. Global stability 

This section is devoted to proving that Reeb graphs of closed curves are stable under 
arbitrary function perturbations. More precisely, it will be shown that arbitrary changes 
in simple Morse functions imply smaller changes in the editing distance between Reeb 
graphs. The proof is by steps: the following Pi'oposition l6.1l shows such a stability property 
when the functions defined on 5' belong to the same arcwise connected component of 
Proposition 16.21 proves the same result in the case that the linear convex combination of 
two simple Morse functions traverses the stratum at most in one point; Theorem |6.3l 
extends the result to two arbitrary functions in 

Proposition 6.1. Let f,g G and let us consider the path h : [0, 1] — > J^(5' ,M) defined 
byh{X) = {l-X)f + lg. Ifh{X) G for every X G [0,1], then d{{rf,f^),{rg,g^)) < 

\\f-8\y- 

Proof. Let 5{h{X)) > be the fixed real number playing the same role of 5 (/) in Theorem 
15.51 after replacing / by h{X). For conciseness, let us denote it by 5{X), and ||/ — g||c2 
by a. If a = 0, the claim trivially follows. If a > 0, let C be the open covering of [0, 1] 

constituted of open intervals 7;^ = — ^^,X + ■ Let C' be a finite minimal (i.e. 

such that, for every /, I^. ^ U ) sub-covering of C, with Xi < X2 < ... < X„ the middle 

points of its intervals. Since C' is minimal, for every i E {1, . . . ,n — I}, I^. f^hi+i is non- 
empty. This implies that 

2 ^ SjXj) 8{Xi+i) max{5{Xi),5{Xi+i)} 

(6.1) Xi+i ~Xi < — ! < . 

Za la a 

Moreover, by the definition of h and the linearity of derivatives, it can be deduced that 

(6.2) \\h{Xi+,)-h{Xi)\\c2 = (A,-+i-A,-)-||/-^llc2- 
Now, substituting (16.1b in ( 16.2b . we obtain 

\\h{h+,)-h{X,)\y < . ||/-g||^, ^max{g(A,),g(A,+,)}. 

Let (r/,(;i^.),/j(Aj)|) be the labelled Reeb graphs associated with {S^ ,h{Xj)), j = 1, . . . ,n. 
Let ; G {l,...,n — 1}. If max{5(A,), 5(A,+ i)} = 5(A,), then using Theorem l5.5l with / 
replaced by /;(A,), g by h{Xi+i) and 5 by ||/!(A,+i) — h{Xi)\\(^, it holds that 

(6.3) fl'((r,(;L,).MA,-)|),(rA(A,^0'M^;+i)|)) < \\h{X,+i)-h{XMc2- 
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The same inequality holds when max{5(/l,), 5(A,+i)} — 5(A,+i), as can be analogously 
checked. 

Now, setting Ao = 0, A„+i = 1, it can be verified that ( 16.31 ) also holds for / = 0,n. 
Consequently, since T/ = ^hCk^}^ and = ^h{X„^^), we have 

«'((r/,j|),(r„^i)) < f rf((r,(;i,),/z(A,)|),(r,(;i,^,)'M^.-+i)|)) < I:iIKA,+i)-KA,)||c2 

i=0 i=0 

= £(A,-+i-A,-)-||/-g||c2 = ||/-^|lc2, 
;=o 

where the first inequality is due to the triangular inequality, the second one to ( 16.31) . the 

n 

first equaUty holds because of ( I6.2l i. the second one because ^ (Aj-|-i — A,) = 1 . □ 

(=0 

Proposition 6.2. Let f^g £ and let us consider the path h : [0, 1] — > ^(5' ,R) defined 
byh{l) = (1 -A)/ + A^. ///!(A) e ^Vo/- every A € [0,1]\{A}, wif/i < A < 1, one//! 
transversely intersects jF' of A, f/;en c/((r^,/j), (rg,g|)) < ||/ — g||c2. 

Proof. We begin proving the following claim. 

Claim. For every 5 > there exist two real numbers A', A" G [0, 1], with A' < A < A", 
suchthatd((r;,(;t')-/<^')|),(rMA")./<^")|)) < 8._ 

To prove this claim, let us first assume that h{X) belongs to To simplify the no- 
tation, we denote /i(A) simply by h. Let p be the sole degenerate critical point for h. It 
is well known that there exists a suitable local coordinate system x around p in which the 
canonical expression of h is h = h{j5) (see Subsection 11.31 and Figure [T] (a) with h 
replaced by /). 

Let us take a smooth function o) : 5^ — > M whose support is contained in the coordinate 
chart around p in which h — h{p) moreover, let us assume that (O is equal to 1 in a 
neighborhood of p, and decreases moving from p. Let us consider the family of smooth 
functions ht obtained by locally modifying h near p as follows: ht = h + 1 ■ (O ■ x. There 
exists f > sufficiently small such that (/) for < t <J, hi has no critical points in the 
support of CO and is equal to h everywhere else (see Figure [T](a) with h, replaced by /2), 
and (//) for — F < f < 0, ht has exactly two critical points in the support of co whose values 
difference tends to vanish as t tends to 0, and ht is equal to h everywhere else (see ||6| and 
Figure[T](a) with /i, replaced by /i). 

Since hi is a universal deformation of h = h{X), and h intersect transversely at 
A, either the maps h(X) with A < A are topologically equivalent to ht with f > or to 
h, with t < (cf. |l6l [m |23]). Analogously for the maps h{X) with A > A. Let us 
assume that h{X) is topologically equivalent to hi with f < when A < A, while h{X) is 
topologically equivalent to /i, with f > when A > A. Hence, for every 5 > 0, there exist 
A', with < A' < A, and A", with A < A" < 1, such that /i(A') and h{X") have the same 
critical points, with the same values, except for two critical points of h{X'), whose values 
difference is smaller than 5, that are non-critical for /i(A"). Therefore, (r/,(;^/),/i(A')|) can 
be transformed into (r/,(;L;/-),/i(A")| ) by an elementary deformation of type (D) whose cost 
is not greater than 5. In the case when h{X) is topologically equivalent to hi with f > 
when A < A, while /z(A) is topologically equivalent to hi with f < when A > A, the claim 
can be proved similarly, applying an elementary deformation of type (B). 
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Let us now prove the claim when h — h{X) belongs to Let us denote by p and 
q the critical points of h such that h{'p) = h{q). Since 'p is non-degenerate there exists 
a suitable local coordinate system x around p in which the canonical expression of h is 
h = h{p) (see Figure [T] (b) with h replaced by /). Let us take (O as before, whose 
support is contained in such a coordinate chart. Let us locally modify h near p as follows: 
hi = h + 1 ■ CO. There exists F > sufficiently small such that for |f| < F, ht has exactly 
the same critical points as h. As for critical values, they are the same as well, apart from 
the value taken at p: hti'p) < h{j}), for — F < f < (see Figure [T] (/?) with ht replaced by 
/i), while ht{p) > h{'p), for < f < F (see Figure[T](/7) with /z, replaced by /2), and ht{'p) 
tends to h{'p) as f tends to (cf. [6^). Since ht is a universal deformation of h — h{X), 
and h intersect transversely at A, we deduce that for every 5 > there exist A', with 
< A' < A and A", with A < A" < 1, such that (r;,(;^/),/!(A')|) can be transformed into 
(r;,(;L//),/z(A")|) by an elementary deformation of type (R) whose cost is not greater than 
5. Therefore the initial claim is proved. 

Let us now estimate d{{r f,f^),{rg,g^)). By the claim, for every 5 > 0, there exist 
0<A'<A"<1 such that, applying the triangular inequality, 

+d{{rj,^^„),h{X")^),{r„g^)) 
< fl'((r^,^),(r,(;^,),/i(A')|))+«'((r,(;L„),/z(A")|),(r„g|)) + 5. 

By Proposition 16. II 

^((r/,^),(r,(,,),MA')|))<||/-KA')|lc2 = A'.||/-g|lc2, 

and 

d((r,(,„),/z(A")|),(r„g|)) < \\h{X")-g\y = (i - a") • Wf-giy. 

Hence, c/((ry,/j), (Fgjg^)) < ||/ — ^||c2 + 5> yielding the conclusion by the arbitrariness 
of 5. □ 

Theorem 6.3 (Global stability). Letf,ge^°. Then d{{T fJ^),{Tg,g^)) < Wf-giy. 

Proof. For every sufficiently small 5 > such that B2{f,5),B2{g,5) C there exist 
fGB2{f,5) andg&B2{g,5) such that the path /i : [0,1] ^{S\R), with h{X) = (1 - 
A)/ + Af, belongs to for every A G [0, 1], except for at most a finite number n of 
values < jUi < ^2 < • • • < M„ < 1 at which h transversely intersects . If n = (« = 1, 
respectively), then the claim immediately follows from Proposition 16. II (Proposition 16.21 
respectively). If n > 1, let < Ai < A2 < . . . < A2„_i < 1, with A2,-i = /X,- for / = !,...,«. 
Then /z(A2,-i) £ for / = 1, . . . ,n, h{X2i) G for /= 1, ...,«- L Set Ao = so that 
/ — h{X()), and A2n = 1 so that ^= h{X2n) (a schematization of this path can be visualized 
in Figure|7|. Then, by Proposition l6.2l we have 

d{{^h{X2,)M^2i)^),{^h{X2,+2)Mhi+2)^)) < \\h{^2i) - h{^2i+2)\\c2 

for every i — Q, . . . ,n — I. Therefore 

di{Tpf^),{T^,g^)) < ";^V(rMA2,),/^(A2/)|),(r,,(;i,,^,),/i(A2,+2)|)) 

< "t\\h{X2d-h{X2i+2)\y < \\f-8\y- 

1=0 
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Then, recalling that / e BiC/, 5) means ||/-/||c2 < S, and B2{f,d) C implies 
that (1— A)/ + A/g for every A G [0,1], we can apply Proposition 16.11 to state that 
d {{T f , f^) , {T p f^)) < d. It is analogous for g and g. Thus, from the triangular inequaUty, 
we have 

d{{rf,f^),{r„g^j) < «'((r^,^),(r^,^))+flf((r^,.^),(r^,^|))+flf((r^,fi),(r„g|)) 
< 28 + \\f-g\y. 

Now, since by the triangular inequality, ||/-§~||c2 < ||/-/|lc2 + Il/-^llc2 + ll^-^llc^' 
with ||/-/||c2<5,and||^-?||c2<5,it follows that «'((r^,j^),(rg,^l))<45+||/-g||c2. 
Finally, because of the arbitrariness of 5, we can let 5 tend to zero and obtain the claim. □ 
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Figure 7 . The linear path used in the proof of Theorem l6.3l 



7. Discussion 

In this paper, we have considered Reeb graphs of curves and have shown that they stably 
represent topological properties of smooth functions. Precisely, we have constructed an 
editing distance between Reeb graphs of closed curves endowed with smooth functions 
/ and g, that is bounded from below by the natural pseudo-distance between {S\f) and 
{S^ ,g), and from above by the C^-norm of f ~ g. 

This paper is meant as a first step toward the study of stability of Reeb graphs of sur- 
faces. While the general technique we use to prove our main result, as well as many 
intermediate results, could be easily generalized to surfaces, the definition of the editing 
distance would need to be appropriately modified. This requires us to classify the possible 
degeneracies of Reeb graphs of surfaces. Moreover, our proof of the metric properties of 
the editing distance exploits some particular properties of curves that are no longer valid 
for surfaces. 

Furthermore, other shape descriptors consisting of graphs constructed out of Morse 
theory, such as the Morse Connection Graph introduced in |9| and further developed in 
lUl, could possibly benefit of some of the results proved in this paper. 

However, some questions remain unanswered also in the case of curves. In the examples 
shown in this paper, the editing distance coincides with the natural pseudo-distance. Is this 
always the case? Moreover, looking at the analogous results proved in ||8] [10| about the 
stability of persistent homology groups, another shape descriptor used both in computer 
vision and computer graphics for shape comparison, we may notice that the C'^-norm rather 
than the C^-norm is used to evaluate function changes. So another open question, strictly 
related to the previous one, is whether it would be possible to improve our result in this 
sense. Other open questions are concerned with applications of the Main Result (Theorem 
16.3b to measure shape dissimilarity coping well with noisy data. On one hand, the result 
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ensures the stability of Reeb graphs against noise, while, on the other, we may wonder 
how likely it is that noise encountered in real data is small with respect to the C^-norm. 
Indeed, it is easy to conceive examples where perturbations that could be seen as noise do 
not correspond to a small value of the C^-norm. For example, the functions represented in 
Figure[8]belong to a sequence of functions (/„) all having the same C^-norm although they 
tend to with respect to the C''-norm. However, one could argue that in a discrete setting, 
at a fixed resolution, sequences of functions as in Figure [8] cannot be found. Moreover, 
this problem would be overcome if the editing distance coincides with the natural pseudo- 
distance. 




Figure 8. The graphs of three functions having the same C^-norm. 
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